Parametric Driving of Dark Solitons in Atomic Bose-Einstein Condensates 
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A dark soliton oscillating in an elongated harmonically-confined atomic Bose-Einstein condensate 
continuously exchanges energy with the sound field. Periodic optical 'paddles' are employed to 
controUably enhance the sound density and transfer energy to the soliton, analogous to parametric 
driving. In the absence of damping, the amplitude of the soliton oscillations can be dramatically 
reduced, whereas with damping, a driven soliton equilibrates as a stable soliton with lower energy, 
thereby extending the soliton lifetime up to the lifetime of the condensate. 

PACS numbers: 03.75.Lm, 05.45.Yv, 42.81. Dp and 47.35. 
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Dark solitons are an important manifestation of 
the intrinsic nonlinearity of a system and arise in di- 
verse systems such as optical fibers '5], waveguides 0, 
surfaces of shallow liquids 4], magnetic films and 
atomic Bose-Einstein Condensates (BECs) 6] . Dark soli- 
tons are known to be dynamically unstable in higher 
than one-dimensional (ID) manifolds (e.g. snake insta- 
bility in 3D systems leading to a decay into vortex rings 
0,11). Solitons in ID geometries experience other in- 
stabilities, whose nature depends on the details of the 
system: For example, dark solitons in optical media are 
prone to nonlinearity-induced changes in the refractive 
index 0] , whereas in harmonically trapped atomic BECs 
they experience dynamical instabilities due to the longi- 
tudinal confinement US d El Eip, [11, as well 
as thermodynamic [Ig and quantum [13 effects. Insta- 
bilities lead to dissipation, which manifests itself in the 
emission of radiation. Compensation against dissipative 
losses by parametric driving has been demonstrated in 
some of the above media il^j . The aim of this Letter 
is to discuss this effect in the context of atomic BECs. 

In atomic gases, the snake instability Q can be sup- 
pressed in elongated, quasi-one-dimensional (quasi-ID) 
geometries and thermal instabilities are minimized 
at very low temperatures T <^Tc (where Tc is the BEC 
transition temperature). In this limit, a dark soliton os- 
cillating in a harmonically confined BEC continuously 
emits radiation (in the form of sound waves) due to the 
inhomogeneous background density _9]. The sound re- 
mains confined and re-interacts with the soliton, lead- 
ing to periodic oscillations of the soliton energy '13| . In 
this Letter we propose the controlled amplification of 
the background sound field, and illustrate the resulting 
transfer of energy into the soliton, in close analogy to 
established parametric driving techniques 0. Energy is 
pumped into the sound field via periodically-modulated 
'paddles', located towards the condensate edge (Fig. 1). 
If the drive frequency is nearly resonant with the soli- 
ton oscillation frequency, one observes significant energy 
transfer to the soliton. In the absence of dissipation, 
this leads to a dramatic reduction in the amplitude of 
the soliton oscillations. Under dissipative conditions, the 



damped soliton equilibrates as a stable soliton with lower 
energy, with its lifetime extended up to the condensate 
lifetime. Moreover, suitable engineering of the phase of 
the driving field (relative to the soliton oscillations) can 
maintain the soliton energy at its initial value for times 
singificantly longer than the undriven soliton lifetime. 

Our analysis is based on the cylindrically-symmetric 
3D Gross-Pitaevskii Equation (GPE) describing the evo- 
lution of the macroscopic order parameter V'(Pi z) of an 
elongated 3D atomic BEC 
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where m is the atomic mass, V — VT{r) + Vu{'r), Vt(i") = 
{m/2){uj'^z'^ + oj'^p^) is the harmonic confining potential 
of longitudinal (transverse) frequency tOz {(^±), where 
^ and VY){r) is the drive potential (Eq. (3)). 
The nonlinearity arises from atomic interactions yielding 
a scattering amplitude g = ATrh?a/m, where a is the s- 
wave scattering length. In this work, a > 0, i.e. effective 
repulsive atomic interactions. The chemical potential is 
given by /i = guQ, where uq is the peak atomic density. 
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FIG. 1: Schematic of parametric driving: (a) Total axial po- 
tential (grey lines, left axis) and density (black lines, right 
axis) of perturbed harmonic trap with propagating dark soli- 
ton (solid lines), at two times (left/right plots) corresponding 
to maximum drive amplitudes, (b) Corresponding drive po- 
tentials (black. Of > 0). Dashed grey lines denote density of 
unperturbed harmonic trap, while dashed black lines indicate 
corresponding density. 
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Dark soliton solutions are supported by the ID form of 
Eq. (1) in the absence of external confinement (V — 0). 
On a uniform background density n, a dark soliton with 
speed V and position (z — vt) has the form, 



ilj{z,t) = V^e^^f^/'*'* i/3tanh 
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(z - vt) 
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where /3 = \J ^ ^ i'^l'^) ' ^'^'^ ^^'^ healing length ^ = 
hj ^fJUm characterises the soliton width. The soliton 
speed v/c = a/1 — (n^/n) = cos(S'/2) depends on the 
total phase slip S across the centre and the soliton depth 
Tid (with respect to the background density), with the 
limiting value set by the Bogoliubov speed of sound 
c — -^//i/m. The energy of the unperturbed dark soli- 
ton of Eq. (2) is given by E^^ = {4:/3)hcn{l - (w/c)2)3/2. 
The drive potential 



Vb = asin(wDi) [e-(^+^°)'/"'o - e-^^"^")'/"-" 
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consists of two periodically modulated gaussian 'paddles', 
with amplitude a, at positions ±zo, oscillating in anti- 
phase at a fixed frequency lou, close to the soliton fre- 
quency (Fig. 1). Such a set-up could be created by time- 
dependent red and blue detuned laser beams with beam 
waist Wo- 

In our work, the dark soliton is defined as the den- 
sity deviation from the unperturbed density (in the 
absence of the soliton) within the 'soliton region' R: 
[zs ~ 5^, Zg + 5^], where Zg the instantaneous position 
of the local density minimum. Such a 'perturbed' 
soliton includes sound excitations located within the 
soliton region at any time. The motion and stabil- 
ity of a quasi-lD dark soliton is well parametrized by 
its energy. In order to facilitate a direct comparison 
of the quasi- ID soliton decay to the analytical homo- 
geneous ID soliton energy, the soliton dynamics are 
parametrized in terms of the 'on-axis' {p = 0) soliton en- 
ergy £'s = {e['!/'(0, z)] - e[V'bg(0, z)]} dz, where £(V') = 

ftV(2m) IWI' + "^^IV'I" + (<7/2) and e[V'bg(0,z)] is 
the corresponding energy contribution of the background 
fluid 13, 14, 20]. 

Dissipationless Regime: Consider first the case of 
no dissipation (Eq. (1)). In the absence of Vd, the soliton 

oscillates at w^oi = w Ja/2 |i,[l3,[il[ll[ll[Il[ll[l3, 

emitting sound waves which oscillate at the trap fre- 
quency uJz ■ This frequency mismatch means that the soli- 
ton propagates through a periodically modulated back- 
ground density, leading to a weak periodic modulation of 
the soliton energy (Fig. 2(a), curves (i)). The amplitude 
of this modulation is enhanced by the coupling between 
longitudinal and transverse degrees of freedom JT, 13]. 

To demonstrate substantial energy transfer into the 
soliton, we start with a low energy shallow soliton (speed 
vq = 0.75c at z = 0). Applying the drive potential in- 
duces an additional periodic background density modu- 
lation and a time-dependence in the soliton oscillation 
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FIG. 2: 'On axis' quasi-lD soliton energy for (i) undriven 
case, (ii) continuous driving and (iii) driving switched off at 
t(j = 80ljJ^, based on simulations of the 3D cylindrically- 
symmetric GPE (black lines) and ID GPE (grey lines) for a 
soliton with initial speed v — 0.75c. (b) Longitudinal soli- 
ton oscillations with continuous driving (grey line), and driv- 
ing switched off at to (vertical grey line) under the 3D GPE. 
Dashed lines indicate corresponding amplitude in absence of 
driving. The trap strength is determined from the chemi- 
cal potential of the system: Quasi-lD: /isd ~ 8to where 
uj = (ujzUj^)^^^ and uj±/ujz = 250. Pure ID: /iiD = IQhuoz 
for which the ID density matches the quasi-lD longitudi- 
nal density. Drive parameters: uu — 0. 98ljso1, Q — O.lfiiD, 
Wo = 3.2lz and zo = 10.7/z where Iz — sjfij (mw-2) the longi- 
tudinal harmonic oscillator length. 



frequency, which is found to vary by no more than 10% 
around its unperturbed value Wsoi- As a result, the rela- 
tive phase between the drive and the soliton oscillations, 
which determines the direction of energy flow between 
soliton and sound, becomes time-dependent. Beginning 
with the drive out-of-phase with the soliton oscillations, 
the soliton initially acquires energy, up to time toi after 
which it begins to lose energy, and the cycle repeats, as 
shown by curves (ii) in Fig. 2(a). This figure illustrates 
the good agreement between the 3D 'on-axis' energy 
(black lines) and the corresponding energy of the pure 
ID simulations (grey lines) (the 3D results feature an ad- 
ditional small amplitude oscillation due to longitudinal- 
transverse coupling). The corresponding beating in the 
soliton oscillation amplitude is shown in Fig. 2(b) (grey 
line) . This beating effect can be visualized as the periodic 
cycling between the initial low energy dark soliton, and 
the nearly stationary high energy soliton. This picture is 
analogous to the cycling of a driven condensate between 
the 'no- vortex' and 'single- vortex' configurations p^ . 

A soliton of higher energy can be created by remov- 
ing the drive potential after a certain pumping time. In 
general, stopping the drive at an arbitrary point in the 
energy gain-loss cycle, will lead to stabilization around 
that energy, with associated residual oscillations. In or- 
der to create a nearly stationary soliton (black line in 
Fig. 2(b)), one must stop the drive when the soliton 
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FIG. 3: Ratio of maximum pumped energy i?f to initial 
soliton energy for a soliton with initial speed uo ~ 0.75c 
in the presence of optimized driving as a function of drive 
frequency (in units of the unperturbed soliton frequency 
Usoi = Wz/\/2) for (i) no damping (black circles), or (ii) with 
damping 7 = 5 x 10~* (grey squares). This value of 7 leads 
to the same soliton lifetime as for the undriven 0.3c soliton in 
Fig. 4. Results based on pure ID GPE for the case of Fig. 2. 



has acquired its maximum energy, for which the energy 
oscillations are suppressed (curve (iii) in Fig. 2(a)). 

The soliton dynamics depend rather sensitively on the 
parameters of the driving field which should be care- 
fully optimized for these effects to be clearly observ- 
able. Firstly, the pumping should take place outside of 
the range of the soliton oscillations (i.e. zq > 9lz for 
vq — 0.75c). If this is not the case, the soliton traverses 
the gaussian bumps, leading to 'dephasing' of emitted 
sound waves and subsequent decay of the soliton [T^. 
The transfer of energy between the soliton and the sound 
field depends on the phase of the drive relative to the 
soliton oscillations, and hence on the drive potential seen 
by the soliton at the extrema of its oscillatory motion. 
This parameter depends in turn not only on the drive 
frequency ljd, but also on the amplitude of the potential 
modulation a, the range of the potential wq, its loca- 
tion zp, and the initial soliton speed vq. If all but one of 
the above parameters are kept constant, then there is a 
resonance around an optimum value of that parameter. 
This resonance is illustrated for the drive frequency wd 
by the open circles in Fig. 3. Note that the maximum 
pumping does not arise at Wgoi, due to the additional fre- 
quency modification induced by the perturbing potential; 
however, the optimum frequency is consistently found to 
lie close to the unperturbed frequency. Importantly, the 
width of the resonance for which the transferred soliton 
energy reaches half its maximum value (FWHM) , is rea- 
sonably broad, of the order of 10% the soliton frequency. 
We have also investigated alternative schemes for pump- 
ing energy into the system. For example, using one off- 
centre paddle, or inverting the sign in Eq. (3) (i.e. a < 0), 
leads to a delayed and less efficient energy transfer, while 
periodically displacing the trap leads to no net increase 
in the soliton energy. 

Dissipative Regime: In a realistic quasi- ID system 
featuring suppressed snake- instability [Sj, both the con- 
densate and the soliton will be prone to damping, e.g. 



due to the presence of a small thermal cloud 0| . A first 
estimate into the effect of dissipation can be obtained by 
introducting a phenomenological damping term H'ydtp/dt 
on the left hand side of Eq. (1). In the absence of driving, 
this term leads to an approximately exponential decay of 
the soliton energy, modified by the oscillatory motion of 
the soliton (dashed lines in Fig. 4(a)). Dissipation also 
damps both the sound field, leading to a narrowing of 
the resonance in the drive frequency, and a shift of the 
resonant frequency towards the unperturbed soliton fre- 
quency (grey data in Fig. 3). 

Stabilization Against Decay: For a high energy soliton, 
continuous parametric driving counterbalances damping 
initially (solid line in Fig. 4(a) up to ojzt ^ 90), but sub- 
sequently the soliton starts to decay (90 < tUzt < 170), 
thereby changing both the amplitude and the phase of 
the soliton oscillations |^, 0| . The evolution in the rela- 
tive phase between the drive and the soliton oscillations 
eventually enables the soliton to gain energy again. Af- 
ter a few such gain-loss cycles, the initial deep i^o = 0.3c 
soliton finally equilibrates as a shallower soliton of energy 
E w 35ftwz (corresponding to vq w 0.7c). 

Stabilization at Fixed Energy: Some applications may 
require a soliton to be maintained at a fixed energy. For 
this, one must balance the competing effects of driving 
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FIG. 4: (a) Soliton energy (initial speed vq = 0.3c) for a dis- 
sipative system (7 = 10""^) in the presence (solid line) or ab- 
sence (dashed) of continuous driving (luo — Usoi, zq = T.llz). 
The driven soliton stabilizes at i5 ~ 35hujz, corresponding to 
Vq ~ 0.7c. (b) Energy of a uq = 0.3c soliton with paramet- 
ric driving (solid black), with a sequence of rephasing cycles 
during the periods between adjacent dashed grey lines. Cor- 
responding curves in the absence of driving are also shown for 
a dissipative(7 = 10"'^, dashed black) and a non-dissipative 
Vq — 0.3c soliton (solid grey), (c) Driven (solid) and undriven 
(dashed) soliton trajectories in the presence of dissipation for 
case (b). Horizontal grey lines indicate the corresponding un- 
damped undriven amplitude of soliton oscillations. Results 
based on the ID GPE, with other parameters as in Fig. 2. 
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and dissipation. This can be achieved by rephasing the 
drive relative to the soUton oscillations at appropriate 
times, as demonstrated by the solid line in Fig. 4(b) 
for a sequence of four rephasing operations. In principle, 
this rephasing can be extended to the duration of the 
condensate lifetime. In an experiment, one actually mea- 
sures the position of the soliton (rather than its energy) 
0. Hence, rephasing could be performed by monitoring 
the amplitude of the soliton oscillations and adjusting 
the drive phase, so that the soliton oscillation amplitude 
remains constant. Fig. 4(c) shows the oscillation am- 
plitude of the parametrically driven soliton (solid black 
line), which is very similar to the undamped undriven 
case (horizontal grey lines), and is clearly distinct from 
the undriven dissipative motion (dashed lines), whose 
amplitude increases until the soliton decays at cuzt = 105. 

Finally, we discuss the relevance of the proposed 
scheme to current experiments with atomic BECs. Given 
a longitudinal confinement u;^ = 27r x 10 Hz, the pre- 
sented results correspond to uj± = 2tt x 2500 Hz and a 
linear 'on-axis' density n = 5 x 10''(1.5 x 10'') m~-^ of ^^Na 
(®^Rb). The harmonic oscillator time unit is r w 15 ms. 
In the dissipative example of Fig. 4 with 7 ~ 10""^, 
this would correspond to a soliton lifetime of around 
1 s, which is consistent with the theoretical predictions 
for solitons in highly elongated three-dimensional ge- 
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H mill- The paddle beams have a waist 
Wo = 20 /im and maximum amplitude a = Ihujz lo- 
cated around zq — 71 z- We have also verified that the 
results presented here hold for smaller aspect ratios (e.g. 
LUj_/u!z — 50). Since, for a given aspect ratio, higher fre- 
quencies correspond to faster timescales, the technique 
presented here is not sensitive to the particular soliton 
lifetime. Increasing, however, 7 beyond an upper limit 
(corresponding to large dissipation), renders it practi- 
cally impossible to pump energy into the system. 

In summary, we have shown that, in the case of a 
dark soliton oscillating in a harmonically trapped Bose- 
Einstein condensate, the addition of two out-of-phase 
gaussian potentials, with amplitude modulated period- 
ically at a fixed frequency close to the soliton frequency, 
pumps energy into the soliton. This technique, which 
bears close analogies to parametric driving, can stabilize 
the soliton against decay for timescales comparable to the 
condensate lifetime. Furthermore, suitable optimization 
of the phase of the drive can force the soliton to maintain 
its initial energy for at least a few times its natural life- 
time. Both effects should be experimentally observable 
in current experiments. 
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